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separated  in  space  and  time,  first  results  show  good  agreement  between  the 
radar,  scintillometer,  and  thermosonde  measurements.^ 
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Sonde  Experiments  for  Comparative 
Measurements  of  Optical  Turbulance 


l.  INTRODUCTION 

Atmospheric  properties  ot  temperatuxe,  pressure,  humidity,  and  -wind  speed, 
as  measured  by  meteorological  radiosondes,  offer  a  means  of  estimating  the  fluc¬ 
tuations  in  the  atmospheric  optical  index  of  refraction.  Fluctuations  of  the  index 

2 

of  refraction,  characterized  as  CM,  are  important  to  numerous  DOD  optical  sys- 

n  2 

terns  being  both  planned  and  developed.  Direct  measurements  of  CN  are  possible 
with  optical  systems  such  as  those  used  to  evaluate  sites  for  astronomical  tele¬ 
scopes.  These  measurements  have  been  limited  to  a  few  locations,  and  this  tech- 

2  2 

nique  is  not  suitable  for  assembling  a  world  wide  data  base  of  C^.  Given  that  CN 
determined  from  meteorological  data  can  be  an  accurate  measure  of  optical  turbu¬ 
lence,  the  best  promise  for  economically  assembling  a  large  data  base  is  through 
the  use  of  the  archived  meteorological  sounding  data  obtained  from  sondes. 

The  present  rawinsonde  measurements  are  usually  achieved  with  altitude 

resolution  of  300  m  or  larger.  This  scale  is  generally  accepted  as  being  larger 

2 

than  tlie  vertical  scales  of  turbulence  that  produce  the  CN’  Thus,  in  order  to  de¬ 
velop  an  analytical  model  predicting  from  the  standard  rawinsonde,  it  is  nec¬ 
essary  to  have  a  statistical  base  relating  the  actual  turbulence  scales  to  the  rawin¬ 
sonde  observed  scales.  New  sondes  (artsondes)  have  been  developed  by  VIZ 
Corporation  for  the  National  Weather  Service  to  measure  the  meteorological 

(Received  for  publication  19  February  1982) 
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variables  at  smaller  altitude  resolution.  These  sondes  were  modified  at  AFGL  to 
2 

provide  data  for  a  Cv  intercomparison  test.  The  program  objective  is  to  compare 
2  n 

as  determined  from  standard  rawinsondes,  the  AFGL  modified  artsondes, 

AFGL  built  thermosondes,  stellar  scintillometers,  aircraft-mounted  commercial 
thermosondes,  and  UHF  radar. 

The  detailed  description  of  the  AFGL  artsonde  and  thermosonde  is  presented  in 
this  report.  Also  presented  is  the  calibration  procedure  and  the  computer  analysis 
used  to  derive  the  meteorological  values  from  the  recorded  telemetry  signal.  The 
sonde  tracking  and  ranging  system  is  described  along  with  the  data  processing  pro¬ 
cedures  used  to  smooth  the  data  and  determine  the  winds  as  a  function  of  altitude. 

o 

A  brief  description  of  experimental  results  is  shown  as  comparisons  of  C^  meas¬ 
urements  obtained  with  radar,  optical  scintillometer,  and  thermosonde. 

2.  RAWINSONDE 

The  rawinsonde  used  in  these  experiments  consists  of  meteorological  measuring 
elements  coupled  to  a  radio  transmitter  and  assembled  Into  a  small  lightweight 

plastic  case.  The  device  is  carried  aloft  by  a  1200-gm  bal'oon  filled  with  helium 
gas.  Included  in  the  train  is  a  small  parachute  to  slow  the  descent  of  the  device 
after  the  balloon  bursts.  An  automatic  launching  reel  is  attached  to  the  train,  which 
lowers  the  radiosonde  until  the  train  is  fully  extended  to  about  2  0  m.  As  the  balloon 
ascends  at  about  5  m/sec,  measurements  of  pressure,  temperature,  and  relative 
humidity  are  transmitted  to  a  ground  station  where  the  data  are  recorded.  A6  the 
balloon  rises,  it  is  followed  electronically  by  radio  direction -finding  equipment  that 
tracks  the  transmitted  signal.  Measurements  of  the  balloon  position  relative  to  the 
ground  station  arc  recorded  in  terms  of  azimuth,  elevation,  and  slant  range.  These 
measurements  are  converted  to  X,  Y,  and  Z  distances.  The  data  are  filtered  as 
described  in  Section  4.4  to  remove  extraneous  signals  and  then  time  differenced 
to  obtain  wind  velocity  profiles. 

Fifty  (50)  state-of-the-art  radiosondes  were  purchased  from  VIZ  Corporation, 
Philadelphia,  Pa.  These  are  called  by  them  the  "Artsonde -Phase  II,  "  and  they 
carry  VIZ  part  number,  1499-311.  This  is  a  1680-MHz  amplitude  modulated  sonde 
with  a  403-MHz  receiver  used  for  trs.nsponder  ranging.  The  transmission  of 
meteorological  (MET)  data  serially  switches  by  solid-state  circuitry  from  pressure 
to  temperature,  to  relative  humidity,  to  reference,  and  back  to  pressure  with  a 
segment  dwell  time  of  2  50  msec.  Premium  sensor  elements  are  used  for  the 
aneroid  baroswitch,  thermistor,  and  hygristor.  Depending  upon  the  sensor  input 
resistance,  a  frequency  of  30  to  1000  Hz  is  produced  that  modulates  the  1G80  MHz 
transmitter.  The  reference  frequency  is  set  close  to  1000  Hz.  It  provides  the 


identification  of  segment  position,  and  it  allows  for  data  correction  of  temperature 
induced  frequency  drifts. 

1  Kawiiisoiuli'  Modifications 

After  delivery  of  the  radiosondes,  sevi-...i  in-house  modifications  were  re¬ 
quired.  Because  greater  temporal  and  spatial  data  resolution  were  desired  than 
were  available  from  a  standard  radiosonde  and  its  associated  GMD  tracker,  special 
payload  and  ground  equipment  changes  were  implemented.  Additionally,  a  sepa¬ 
rate  instrument  to  measure  small  temperature  fluctuations  (thermosonde)  was 
attached  physically  and  electrically  to  the  radiosonde.  The  modification  included 
changes  in  the  radiosonde  transmitter,  receiver,  solid-state  switching,  and  MET 
oscillator  circuit,  boards. 

2.1.1  MET  CIRCUIT  CARD 

In  order  to  transmit  the  additional  thermosonde  data,  the  number  of  commuta¬ 
tion  data  segments  was  Increased  from  the  standard  four  to  eight.  To  accomplish 

this,  a  second  quad-bilateral  switch  (p/n  406(1)  was  added  to  the  MET  circuit 
board.  This  was  facilitated  oy  the  existing  4017  If  clock,  which  allows  serial 
switching  of  up  to  ten  segments.  Two  of  the  additional  segments  were  used  for 
thermosonde  information  and  two  were  used  as  auxiiliary  segments.  The  auxiliary 
segments  allowed  super -com mutation  uf  the  MET  data  or  the  reporting  of  electron¬ 
ics  temperature. 

To  accomplish  calibration  and  ground-base  checks  the  eight  lock-up  and  one 
common  test  points  were  brought  to  an  eight-position  rotary  switch.  During  cali¬ 
bration  the  normally  free-running  segment  switch  was  stopped  on  a  segment  cor¬ 
responding  to  the  position  of  the  switch.  The  frequency  of  each  segment  was 
examined,  in  turn,  by  rotating  the  switch  from  one  position  to  the  next.  The 
measured  frequency  was  then  recorded  along  with  the  measured  temperature, 
humidity,  pressure  commutator  position,  and  thermosonde  test  value.  The  baro- 
switch  that  forms  one  part  of  the  MET  board  transmits  one  of  four  unique  fre¬ 
quencies  depending  upon  the  position  of  the  baroswitch  armature.  The  original 
separation  of  frequencies  was  only  10  or  20  Hz.  In  order  to  provide  a  separation 
of  100  Hz  between  each  frequency,  three  new  resistors  replaced  existing  ones  with 
appropriate  values.  This  reduced  confusion  over  contact  identification. 

2.1.2  TRANSPONDER  RECEIVER 

The  balloon  borne  403-MHz  transponder  receiver  is  used  to  derive  the  slant 
range  in  the  following  manner:  Ground-base  equipment  generates  a  PCM  code  that 
is  synchronized  with  the  start  pulse  of  a  time  interval  counter.  This  code  modulates 
the  <i03-MHz  uplink  transmitter.  The  receiver  transponder  uses  this  code  to 


modulate  a  70-kHz  subcarrier  oscillator.  This  subcarrier  is  mixed  with  the  MET 
data  and  used  to  modulate  the  balloon-borne  transmitter.  Ground-base  equipment 
retrieves  and  decodes  the  PCM  code,  which  then  sends  a  stop  pulse  to  the  inter¬ 
val  counter.  Since  slant  range  is  calibrated  vs  interval,  a  direct  digital  coded 
range  output  is  produced. 

2.  l.  3  TRANSMITTER 

Normally  the  transmitter  is  100  percent  amplitude  modulated  bv  .ne  MET  data 
pulses.  When  the  MET  data  pulse  is  high,  the  transmitter  Is  'hut  down.  This 
results  in  an  AM  pulse  repetition  frequency  that  is  measured  by  ground  receiving 
equipment  to  determine  the  sensor  input  resistance. 

The  on-board  403-MHz  receiver  frequency  modulates  the  transmitter.  Since 
this  information  is  lost  when  the  MET  signal  is  high  (because  the  transmitter  is 
shut  down),  difficulty  in  holding  synchronization  to  the  ranging  signal  develops. 
Consequently  the  evaluation  of  slant  range  data  contains  associated  errors.  To 
eliminate  this  problem,  the  ac-coupllng  AM  input  to  the  transmitter  was  removed. 
Instead,  the  MET  data  was  ac-coupled  to  the  transmitter  In  the  same  manner  and 
to  the  same  junction  point  as  the  ranging  signal.  This  results  in  MET  data  fre¬ 
quency  modulation  of  30  to  1000  Hz  and  ranging  modulation  of  70  kHz. 

2.2  ReauUnce  Calculation 

This  section  presents  the  MET  data  transfer  equations  as  given  by  VIZ 
Corporation  (document  nos.  VIZ  800721  and  761101B).  The  form  of  these  ex¬ 
pressions  as  well  as  the  constants  depend  on  the  sensor  type  and  manufacturing 
lot  number.  The  input  resistance  to  the  radiosonde  corresponding  to  the  measured 
MET  frequency  is: 

R input  Rlnt  <fref/fmeas  "  ' 


where  f  ^  is  the  average  measured  reference  frequency  and  fme&s  is  Hie  average 
measured  segment  frequency  when  the  baroswitch  in  on  a  space  (that  is,  between 
two  discrete  contact  points). 

The  thermistor  or  pressure  resistance  is  the  input  resistance  when  on  a  tem¬ 
perature  or  pressure  segment. 

Tlie  input  resistance  when  on  a  humidity  segment  is  the  equivalent  parallel 
resistance  of  a  1-M  resistor  and  the  hygristor  resistance;  thus. 
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2.2.1  TEMPERATURE  CALIBRATION 

Temperature-  is  calculated  from  e.c  thermistor  resistance  K,p  as  follows: 


-  273.  15 


i  Ak  i'u  niT/R30n! 


where  R„,,  ts  the  acculok  lock-in  resistance  at  30°C.  The  coefficients  for  preini. 
ov 

temperature  element  lot  PG97n  as  provided  by  VIZ  are: 


A  =  3.2987  E-03 
o 

A  =  4. 7610  E-04 
A2  =  2.8417  E-OS 
A3  =  1.  5691  E-06  . 


2.2.2  HUMIDITY  CALIBRATION 


Relative  humidity  is  calculated  from  the  hygristor  resislan  oe  RRH  as  follows 


RH  -  A  - 


I  Dk«f(RRH>lk 


where 


f^KRH)  =  ln  ^?RH^R33^  ' 

R33  is  the  acculok  lock-in  resistance  at  33  percent  RH  and  2S°C 


fit)  =  Y  cK  (DV 


(t  =  °Celsius)  . 


The  coefficients  fo"  premium  humidity  element  lot  22G5  as  provided  by  VIZ  are: 


for  *WR33  ~1: 

C  =  7.885  E-01 
o 

C  =  9.286  E-03 
C2  =  -2.  462  E-05 
C3  =  -3. 368  E-07. 


for  rhh/r33  <  1: 

C  =  9.243  E-01 
o 

=  3. 059  E-03 
C  =  -1.  188  E-06 
C3  =  0. 
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for  f(RRH)  2:  -0.2: 

Do  “  l- 

Dj  «  7.  1454  E-01 
D2  *  -7.9163E-02 
Dg  «  3.2555  E-02 
D4  «  3.9255  E-03 
A  «  102. 


for  f(RHH)  <  -0.2: 

D  «  -2.09Q9E00 
0 

Dj  ■  9.2406  E00 
D2  =  2.  1693  E01 
Dg  =  2.2214  E01 
D4  *  -6.  1926  E00 
A  ■  0. 


2.2.3  PRESSURE  CA LIBRA TION 


Each  radiosonde  is  provided  with  pressure  calibration  values  for  the  leading 
edge  of  160  contact  positions  of  the  baroswitch  by  VIZ  Corporation.  Four  discrete 
frequencies  can  be  transmitted  by  the  baroswitch,  depending  upon  at  which  of  the  four 
classifications  of  contacts  the  baroswitch  is  located.  The  four  frequencies  were 
adjusted  to  provide  approximately  500“,  600-,  700-,  and  800-Hz  signals.  Approxi¬ 
mately  500  Hz  is  transmitted  when  the  baroswitch  is  located  between  any  two  dis¬ 
crete  contact  points.  Approximately  600  Hz  is  transmitted  for  the  following  contact 
numbers:  30,  45,  60,  75,  90,  105,  120,  135,  140,  145,  150,  155,  and  160.  These 
start  at  contact  number  30  and  increase  by  15  for  contact  numbers  less  than  135  and 
increase  by  five  foi  contact  numbers  greater  than  135.  Approximately  700  Hz  is 
transmitted  for  the  following  contact  numbers:  5,  10,  15,  20,  25,  35,  40,  50,  55, 
65,  70.  80,  85,  95,  100,  110,  115,  125,  130,  136,  137,  138,  139,  141,  142,  143, 
144,  146,  147,  148,  149,  151,  152,  153,  154,  156,  157,  158,  and  1.50.  These  start 
at  contact  number  5  and  increase  by  five  for  contact  numbers  less  than  135  and 
increase  by  one  for  contact  numbers  gi  eater  than  135.  Approximately  800  Hz  is 
transmitted  for  all  other  contacts.  Before  each  launch  the  baroswitch  is  manually 
set  to  the  contact  position  corresponding  to  the  ambient  station  pressure.  During 
analysis  the  computer  identifies  groups  and  subgroups  of  contacts  and  defines  a 
particular  contact  within  a  larger  group.  Thus,  if  transmission  is  lost  for  a  short 
time,  the  computer  redefines  the  location  of  the  baroswitch  and  estimates  the  pres¬ 
ent  contact  number.  Pressure  is  defined  at  the  leading  edge  of  each  contact  and  is 
logarithmically  interpolated  between  contacts  based  on  altitude. 

Sample  results  of  the  pressure,  temperature,  humidity,  and  calculations 
are  presented  in  Figure  1.  The  first  column  is  the  record  number.  "Time"  is 
defined  as  seconds  after  launch.  "C'TKMP"  is  the  corrected  temperature  based  on 
inferred  time  lags.  Flo  correction  is  made  to  humidity  below  -40°C.  Pressure 
contacts  are  split  into  classes  0  to  3.  Column  11  shows  that  a  contact  number  is 
defined  at  the  first  encounter  of  a  new  contact.  Columns  14  and  16  show  that 
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Figure  1.  Sample  CXitput  of  Heduced  Kadio^onde/Fhormosonde  Data.  Processed  results 


(AT)rms  is  undefined  for  values  less  than  0.002°.  The  "  in  the  last  column 

flag  values  of  (AT)rmg  low  gain  that  are  greater  than  0.  035°.  The  corresponding 

high  gain  value  is  not  valid  for  this  condition. 


3.  THERMOSONDE 


The  thermosonde  is  a  small,  lightweight,  temperature-sensing  electronic 

unit  that  attaches  to  the  radiosonde.  It  measures  the  difference  In  temperature  of 

the  ambient  atmosphere  between  two  probes  mounted  horizontally  1  m  apart.  As 

ihe  balloon  ascends  or  descends,  the  temperature  differences  are  processed  and 

routed  to  the  radiosonde  MET  segments.  The  received  data  is  used  to  calculate 

2 

the  atmospheric  temperature  structure  parameter  C^,.  These  data,  together  with 

pressure  and  temperature,  are  used  to  calculate  the  refractive  index  structure 

2  2  1 
parameter  C^.  Fundamentally,  is  defined  as: 

C^j  *  [n(f,)  -  n(r“)]1 2 3/r2/3  . 


for  n(rj)  the  index  of  refraction  at  point  rt  and  r  =  |  rj  -  r2  I  •  This  parameter  is  the 
key  unknown  in  estimating  the  performance  cf  laser  systems  through  the  clear 
atmosphere.  Among  thoce  explicit  or  indirect  dependencies  on  the  profile  are: 
log-amplitude  of  light  fluctuations,  wave  structure  function,  receiver  diameter  for 
heterodyne  detection,  and  modulation  transfer  function. 

In  flight,  the  thermosonde  passes  its  temperature  difference  measurement  of 
T2  -  T^  through  a  root  mean  square  (rms)  module.  Ground-based  process!'  g  con¬ 
verts  this  to  the  mean  square  temperature  fluctuation  (T,,  -  T~P ,  which  over  one 

2  •)  &  1 

meter  separation  equals  This  is  converted  to  C^,  through  the  equation: 


C^(h)  = 


79.  9  P(h)  X  10 


-C 


T2(hi 


x  c: 


(1) 


where  P(h)  is  pressures  in  mbar  and  h  is  altitude.  * 

The  original  lightweight  thermosonde  system  was  developed  by  GTE  Sylvanla 

2 

lr.c.  ,  Electro-optics  Organization  for  the  Goddard  Space  Flight  Center  in  197  ].  ' 

1.  Tatarski,  V.I.  (196 1)  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill 

Co.  ,  New  York. 

2.  'Iitterton,  P.J.,  Mallery,  L.E.,  and  Arken,  T.A.  (197 1)  Eight  weight 

Thermosonde  System,  GTE  Sylvanla,  Inc.,  Final  licport,  NATiJT- 1 14 93 . 

3.  GTE  Sylvanla,  Inc.  (1972)  Lightweight  The rmosonde  System  Model  HO 

(Instruction  Manual),  Mountain  View,  California,  IJ440. 
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figure  2,  Photograph  of  Thermosonde  and  Micro  Thermal  Probe.  Tungsten  wire 
diameter  3.45  *im.  Effective  length  Is  0.45  cm.  Nominal  ice  point  resistance  ie 
27  ohms  and  its  temperature  coefficient  of  resistance  is  0.  00375°C,'_1.  , Sensor  wire 
retouched  for  better  photographing 


3000  II z  and  0.2  llz,  respectively.  A  1-sec  rms  average  ie  performed  on  the 
filtered  signal. 

In  the  ease  of  a  square  wave  signal,  the  rms  reeistance  ie  one-half  the  peak- 
to-pcak  resistance  and  Eq.  (2)  becomes 

AT  »  AH  / (2  at  R  )  ,  (3) 

rms  p  o 

The  fine  V'J.  08  percent  pure  tungsten  wire  in  manufactured  by  the  Sigmond 
Cohen  Corp.  ,  Mt.  Vernon,  N.Y.  ,  and  the  wire  iB  attached  to  a  probe  assembly 
by  Outunietrlcs,  Inc.  The  nonil  .al  wire  diameter  is  3.  45  pm.  A  0.45-cm  length 
of  the  wire  is  mounted  between  two  0.  5-ern  long  wire  supports  extending  out  from 
the  probe  holder  (f  igure  2).  Its  nominal  ice  point  resistance,  I?o,  is  27  ohms  at 
Cr’c.  Thin  sets  the  electrical  resistivity  at  5.0  nil  -  cm.  'The  temperature  coef¬ 
ficient  of  resh  tance  wan  measured  at  0.  0037 U°V  .  Note  that  the  electrical 


Pi 


HTTP*. 


resistivity  and  temperature  coefficient  differ  from  tabulated  values.  This  is  pro- 

8 

bably  due  to  a  low  annealing  temperature  of  the  drawn  wire. 


U  Swot  H»iHi|  Effects 

Optimum  performance  of  the  thermosonde  can  be  achieved  by  passing  a  current 
through  the  wire  that  is  large  enough  to  enhance  the  temperature  response,  but 
small  enough  to  avoid  significant  velocity  response.  Since  a  cool  ambient  fluid 
flowing  over  a  warmer  body  will  lower  the  temperature  of  the  body  toward  the 
fluid  temperature.  It  Is  desirable  to  keep  the  current-induced  temperature  dif¬ 
ference  below  the  minimum  detectable  value  of  0.  002°C  rms.  The  maximum  cur¬ 
rent  allowed  for  not  exceeding  this  temperature  difference  can  be  estimated  from 
a  heat  transfer  analysis.  The  current  generated  temperature  difference.  AT^. 
can  be  determined  as 


ATW 


OP 
JA  h 

S 


(4) 


where 

OP  is  the  electrical  power  input  between  two  conditions 

J  is  the  mechanical  equivalent  of  heat  4.  1.8  (Joule/Cal) 

-4  2 

A„  is  the  surface  area  of  the  wire  =  5  X  10  (cm  ) 

“  _tj  ___  A 

h  is  the  heat  transfer  coefficient  (Cal  *  s  *  cm  *  K  ). 

The  power  input  is  P  =  I^(RA  +  OR)  where  RA  is  the  sensor  resistance  in  the  ab¬ 
sence  of  an  applied  current  and  OR  is  the  change  in  resistance  due  to  applied 
current. 

The  heat  transfer  coefficient  of  the  actual  thermosonde  probe  was  determined 
in  the  laboratory  by  the  temperature  difference  between  two  power  levels.  A  spe¬ 
cial  bridge  te3t  circuit  (Figure  3)  was  built  to  provide  a  constant  current  across 
the  bridge.  This  current  alternatively  can  be  stepped  high  and  low,  depending  upon 
the  input  voltage  E.  The  bridge  output,  E^,  was  amplified  by  a  factor  of  83.  33. 

The  additional  wire  resistance,  R,  caused  by  the  current  heating  was  determined 
from  the  observed  voltages  and  current  as 

11  =  Eo/83.33Ij  .  (5) 


The  current  converter  input  voltages  were  stepped  from  0.  5  V  to  5  V,  which  pro¬ 
duced  probe  currents,  L,  of  0.  5  mA  and  5  rnA,  respectively.  The  observed 


8.  SmitboMs,  C.  J.  (1952)  Tungsten,  A  Treatise  on  Its  Metallurf 
And  Applications,  Chapman  &  Hall,  London,  3rd  ed. 
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Figure  3.  Schematic  of  Thermosonde  Sensor  Convective  Heat  Transfer  Test 

Circuit.  For  a  given  input  voltage,  E^n,  the  LH0041  circuit  produces  a  con¬ 
stant  probe  current  of  (Ein/(2  X500))  amps.  The  725  amplifier  provides  an  out¬ 
put,  Ee  ■  83.33  X  Bridge  potential  difference.  Thus  the  differential  resistance 
of  the  bridge  due  to  probe  I^R  heating  is:  R  «=  EQ/83.  33  Ij, 


outputs  were  0.  01  V  and  0.  32  V.  The  temperature  difference  produced  in  the 
wire,  is 

ATW  =  ^5.  0  '  R0.  5^  “  Ro  *  (6) 


where  in  this  case  R^  =  28.28  n,  AR  =  Rg  Q  -  R0  g  =  0.  528  fi,  and  AT^,  ■  5.  36°C. 

The  difference  in  power  between  the  two  test  conditions  was  determined  as 
AP  ■  P,.  q  -  Pq  5  =  (24.75  R^  +  60  Eg  Q-  6EQ  g)  pW  with  a  measured  RA  =  28.25  SI, 
AP  *  718  pW,  From  Eq.  (4),  the  experimentally  determined  heat  transfer  coef¬ 
ficient  is  found,  h  =  0.  064  cal  •  s  *  •  cm  2  •  °K 

In  flight,  using  a  probe  power  of  l / 3  pW  rms  produces  a  self-heating  temper¬ 
ature  rise  [using  Eq.  (4)]  of 

AT  =  0.  002 5°C  (rms)  . 


Consequently,  self -heating  effects  set  the  thermosonde  noise  level  at  0.  002°C.  The 
actual  flight,  condition  heat  transfer  coefficient  depends  upon  air  velocity  and  den¬ 
sity.  A  study  was  performed  to  determine  the  convective  heat  transfer  coefficient. 
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as  a  function  Reynolds  number.  The  heat  transfer  coefficient  h  is  described  in 
terms  of  the  Nusselt  number  Nu^  as  h  =  Nu^k/d,  where  k  is  the  thermal  conductivity 
of  the  air  and  d  is  the  characteristic  length.  The  Nusselt  number  can  be  calculated 

q 

from  the  Reynolds  number  Re^  over  limited  regions.  Calibrations  of  h  vs  pressure 
and  wind  velocity  were  performed  in  a  bell  jar/blower  vacuum  system  at  conditions 
corresponding  to  the  surface  level,  that  is.  p  =  760  nimHg,  T  =  24°C, 
velocity  =  5m  •  s  1  (corresponding  to  the  vertical  balloon  ascent  rate)  Re  =  1.2; 
h  is  evaluated  to  be  0.  12  cal  •  s  *  *  cm  2  •  °K  *  or  double  the  "no-wind”  condition. 
At  lower  pressures,  the  bell  jar  density  is  correlated  with  the  standard  atmosphere 
to  obtain  altitude.  An  altitude  profile  of  h  is  constructed  from  these  calibrations 
(Figure  4).  Above  29-km  altitude  molecular  flow  occurs.  Here,  the  free  molecular 
heat  conductivity  given  by  Dushman  is  used.^  It  is  clear  that  an  increase  in  velo¬ 
city  (Reynolds  number)  at  standard  conditions  enhances  the  value  of  h  and  lowers 
self-heating  even  below  the  specified  0.  002°C.  At  the  maximum  expected  altitude 
of  about  30  km,  h  *  0.031  cal  •  s"1  •  cm"2  *  “K"1,  and  self-heating  can  be  ex¬ 
pected  to  approach  a  noise  value  of  0.  005°C. 


Figure  4.  Heat  Transfer  Coefficient  Curve  of  Therm o- 
sonde  Sensor  as  a  Function  of  Altitude.  Bell  jar  cali¬ 
bration.  At  STP,  h  is  measured  at 
0.064,  cal  •  s’1  •  cm"2  •  °K'l.  WUh  air  speed  = 

5  m  ‘  8"  h  =  0.  12  oe.1  •  s’l  •  cm"*-  °K'1  at  30  km, 
h  =  0.  031  cal  •  b“*  ■  cm~2  •  Above  30  km, 

free  molecular  conductivity  applies 


9.  Marton,  L.  (1981)  Fluid  Dynamics.  Academic  Press,  New  York,  Vol.  18, 

Part  A,  p.  271. 

10.  Dushman,  S.  (1962)  Scientific  Foundations  of  Vacuum  Techniques.  John  Wiley, 
New  York,  p.  47.  2nd  ed. 
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3.3  Senior  Frequency  Reeporue 

If  the  temperature  of  the  environment  suddenly  changes,  one  wishes  to  know 
how  fast  the  sensor  will  follow  the  environmental  changes.  Assume  that  the  wire 
Initially  has  a  temperature  and  Is  to  be  brought  to  an  environment  temperature 

T  .  The  time  for  the  wire  temperature  to  recover  a  given  fraction  of  the  environ- 

e  11 
mental  temperature  Is  given  as: 


%  =  .1 .2  abt. 

V  '  r 


(7) 


where, 

% 

6. 


r 

a 

b 

h 

k 


pw 


a 

b 

t 


max  e 

VTe 

wire  radius  (cm) 

n 

thermal  dlffusivlty  of  die  wire  (cm  /sec) 
h/k 


-1 


,-2 


°K"1) 


the  maximum  temperature  in  the  wire  after  time,  t 


thermal  conductivity  of  tungsten  {cal  •  6 


-1 


r1  •  °k-1) 


-i 


-i 


°K-i 


0.426  cal  •  s  "  •  cm 
k/<'wCpw) 

tungsten  density  *>  19.  3  g  '  cm 
heat  capacity  of  tungsten  =  0.  032  cal  •  g 


-1 


V1 


0.  69  cm  ’  8 
•1 


-1 


!  0. 16  cm 
!  time  (sec) 

Thus,  the  response  time  required  to  recover  0o/<h  "  percent  and  =  90  per¬ 

cent  is  0.4  msec  and  0.72  msec,  respectively. 

This  was  verified  experimentally.  The  sensor  was  placed  in  a  bell  jar  vacuum 
chamber.  A  current  pulse  train  of  0.  5-  to  5-  mA  amplitudes  was  applied  to  the 
sensor  and  the  ac  bridge  signal  response  was  displayed  on  an  oscilloscope.  It  was 
observed  that  the  bridge  output  approached  l/2  of  its  full-scale  change  in  0.  3  msec. 
This  gives  a  frequency  response  f ■  3570  Hz.  It  Is  clear  that  the  sensor  has  an 
adequate  temperature  frequency  response  for  the  purpose  of  turbulence 
measurements. 


11.  Jakob,  M.  (1949)  Heat  Transfer.  John  Wiley,  New  York,  Vol.  1,  p.  291. 


20 


3.4  Thermosonde  Calibration 

The  thermosonde  measurement  is  transmitted  as  a  frequency  value,  nominally 
ranging  between  30  and  700  Hz,  Each  frequency  is  related  to  a  resistance  which 
is  defined  by  the  observed  temperature  difference  using  Eq.  (3).  The  relationship 
between  frequency  and  resistance  is  determined  by  actual  observations.  A  thermo - 
sonde  test  box  was  fabricated  that  alternately  switches  a  test  resistor  in  and  out 
across  one  of  the  pair  of  matched  27  fl  resistors  to  simulate  the  probes.  Any  one  of 
12  precision  resistors  ranging  from  1.4  M  n  to  20  k  O  can  be  inserted  for  the  test 
resistor.  These  parallel  resistors  establish  a  known  equivalent  resistance.  The 
differential  resistance  (AR)  across  the  bridge  is  then. 


AR  «  R  - R 

ps  eq 


where  Rpg  is  the  simulated  probe  resistor  *  27  Cl,  and 


R  XR 

r,  _  D8  8 

Req  R  “  +TT 

ps  s 


where  R  is  the  inserted  test  resistor.  Thus, 
s 


AR 


£®_ 
+  R 


s 


For  each  thermosonde  a  table  is  constructed  for  AR  vs  radiosonde  frequency. 

ATrmig  i3  calculated  from  Eq,  (3)  and  a  polynomial  fit  is  performed  on  this  data 
in  order  to  provide  the  transfer  function.  Typical  calibrations  are  shown  in  Fig¬ 
ures  5  and  6. 

3.5  Thermownde  Environment 

The  thermosonde  was  suspended  18  m  below  the  sonde  balloon.  The  balloon 
diameter  expands  from  approximately  2  m  at  launch  to  about  10  m  at  30  km. 

A  working  rule  is  to  locate  experiments  at  least  five  diameters  below  the  bal¬ 
loon  to  avoid  wake  effects.  The  18 -m  suspension  length  is  a  compromise.  A 
longer  length  can  enhance  a  pendulum  motion  that  would  interfere  with  measure¬ 
ments.  Thus,  it  is  expected  that  the  thermosonde  will  be  located  inside  the  balloon 
wake  except  during  instances  of  large  shears.  However,  even  with  the  thermosonde 
inside  the  balloon  wake,  it  is  not  clear  what  type  of  thermosonde  error,  if  any,  the 
wake  would  produce. 
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Figure  5.  Calibration  Curve  of  Thermosonde  Sensor 
AR  as  a  Function  of  Radiosonde  Frequency.  Test 
resistors  are  switched  quickly  across  simulated 
probes.  This  provides  a  calibration  frequency  for 
each  known  differential  resistance 


FREQUENCY (Hii 


Figure  6.  Calibration  Curve  of  Thermosonde 
(AT)rrns  as  a  Function  of  Radiosonde  Frequency.  For 
a  probe  of  given  resistance  and  temperature  coeffi¬ 
cient,  the  resistance  vs  frequency  calibration  curve 
is  converted  to  (AT)rmB  vs  frequency 
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The  existence  of  a  wake  behind  a  slow  moving  balloon  (Re  ~  10  )  could  gener¬ 
ate  a  temperature  fluctuation  due  tc  two  mechanisms.  The  first  is  due  to  turbulent 
mixing  of  the  atmospheric  temperature  structure  over  some  vertical  scale  length. 
The  second  could  arise  as  a  companion  to  compressible  flow  fluctuations  of  the 
velocity  in  the  wake. 

In  the  troposphere,  the  sonde  is  usually  in  supercritical  flow,  and  the  sepa¬ 
rated  flow  produces  organized  vortices.  At  about  11-kms,  the  flow  becomes  sub- 
critical  and  the  flow  around  the  balloon  turns  into  a  randomized  separation  and  a 
uniform  turbulent  wake.  Analysis  of  the  thermosonde  data  reveals  no  significant 

differences  between  tropospheric  and  stratospheric  data.  The  thermosonde  data 

2 

indicates  the  probability  of  C_  is  log  normally  distributed.  Power  spectral  analy- 
2  1 

sis  of  C^,  shows  that  no  frequency  Spikes  or  features  exist  over  the  measurement 

range  of  2-  to  200-sec  periods.  This  implies  that  any  wake  effects  arising  from 

atmospheric  structure  would  have  to  arise  from  scales  smaller  than  10  m.  The 

2 

thermosonde  observes  layers  of  high  CT  over  scales  larger  than  2  0  m  as  well  as 
2  ^ 

undetectable  C'X,  (S/N  <  1)  for  scales  larger  than  10  m.  This  strongly  suggests 
that  the  existence  of  a  wake  does  not  cause  false  signals. 

The  explanation  of  why  there  is  no  wake  effect  could  be  due  to  the  rapid  dis¬ 
sipation  of  wake  gradients.  The  small  scale  of  the  wake  (length  less  than  10  m) 
allows  rapid  mechanical  mixing  to  a  still  smaller  scale  followed  by  rapid  thermal 
dissipation.  We  have  made  measurements  of  the  rate  of  energy  dissipation,  e,  in 
the  wake  of  high  altitude  aircraft.  Film  density  analysis  of  aircraft  wakes  marked 
by  creating  a  smoke  trail  of  titanium  oxide  hydrates  shows  the  energy  decay  rate 
to  follow  at4  dependence,  Thus,  for  a  vertical  rise  rate  of  5  m/sec,  the  balloon 
wake  energy  should  be  expected  to  dissipate  to  1/256  of  its  original  value  by  the 
time  it  reaches  the  thermosonde. 

Confidence  in  the  absence  of  wake  effect  will  be  developed  as  more  of  the  data 
2 

is  analyzed  and  comparisons  are  made  with  the  scintillometer  and  radar. 

4.  DATA  REDUCTION 

The  "raw"  transmitted  MET  signal  is  bandpass  filtered  at  400  Hz  and  5000  Hz, 
after  which  it  triggers  a  "one  shot"  comparator.  This  results  in  a  "clean"  train 
of  constant-amplitude/constant-width  pulses  equivalent  to  the  MET  frequencies. 

A  Hewlett-Packard  model  5226b  timer/counter  measures  the  time  period  of  every 
other  frequency  period  and  transfers  these  times  in  a  binary  coded  decimal  format 
to  a  Honeywell  computer,  model  3 16.  The  computer  packs  a  magnetic  tape  eight 
bits  per  word  in  a  pre -determined  record  size  and  format.  These  tapes  are  then 
unpacked  and  converted  to  CDC  6600  computer  compatible  9  track  tapes. 
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4.1  Tine  Period  Record* 


^  Running  time  ie  determined  by  setting  the  algorithm  time  =  Header  time  + 

2  i  T  ,  where  header  time  is  the  actual  Zulu  time  at  the  start  of  each  compu- 
1« !  true 

ter  record.  The  factor  2  appears  because  the  counter  measures  every  other  peri¬ 
od.  The  true  time  period  is: 


T 

true 


tape  speed  input  to  counter 
tape  speed  of  original  recording 


X  measured  counter  time  . 


For  this  series,  the  tape  speed  counter  input  was  l/2  the  original  tape  speed.  If 
any  T  true  is  less  than  0.  91  msec,  it  is  considered  a  data  glitch  and  dropped  from 
the  summation  and  analysis. 


4.2  Reference  Frequency  ID 

The  unpack  routine  is  initialized  to  the  measured  reference  frequency  at  launch. 
A  forward  test  searches  for  three  successive  data  periods  with  a  frequency  of  the 
Initialized  value  ±  5  percent.  When  this  criterion  is  achieved,  a  backwards  test 
searches  in  reverse  from  the  first  of  the  three  found  frequencies  to  the  first  en¬ 
counter  of  the  initialized  reference  frequency  ±50  Hz.  Once  aloft,  the  reference 
frequency  is  allowed  to  drift  slow1”  *ora  its  initialized  value. 


4.3  Segment  Window  ID 

After  locating  the  last  reference  period  of  the  previous  reference  segment,  the 
time  interval  of  the  entire  eight  segment  commutation  frame  is  measured  by  sum¬ 
ming  each  period  and  multiplying  the  result  by  two  as  described  before.  This  de¬ 
fines  the  frame  time.  Since  there  are  eight  segments,  the  frame  time  is  divided  by 
eight  to  identify  each  segment  interval.  Since  the  frame  rate  is  set  at  approximate¬ 
ly  2 .  0  sec,  the  segment  dwell  time  is  approximately  2  50  msec  per  segment.  Only 
the  center  60  percent  of  each  time  segment  is  used  in  the  frequency  determination 
for  that  segment.  The  first  and  last  20  percent  of  the  segment  time  is  dropped  and 
the  remaining  periods  are  averaged.  Inverting  the  average  period  provides  the 
average  frequency  for  the  segment.  The  standard  deviation  of  each  segment  inter¬ 
val  is  calculated,  and  if  a  -  exceeds  ±  10  Hz,  or  if  cn  __  ,  a.,u,  or  ct  exceeds 

rel  temp  r(H  press 

±  5  Hz,  the  two  smallest  time  periods  are  dropped.  In  most  cases  this  removes 
spurious  pulses,  A  new  average  and  o  are  redetermined,  and  if  aref  exceeds 
±  10  Hz,  the  entire  commutation  frame  is  dropped  from  analysis.  If  a^enip>  cthh’ 
or  °press  exceeds  ±  5  Hz,  the  calculated  segment  frequency  is  considered  erron¬ 
eous,  and  individual  bad  segments  are  dropped  from  analysis.  In  addition,  if  at 
least  two  time  periods  are  not  found  in  a  particular  segment,  that  segment  is 
dropped  from  analysis.  The  segment  sequence  per  frame  is,  in  most  cases: 
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reference,  thermosonde  high  gain,  temperature,  relative  humidity,  thermosonde 
low  gain,  pressure,  temperature,  and  relative  humidity.  When  the  thermosonde  is 
not  flown,  those  segments  usually  contain  second  samples  of  pressure  and  relative 

humidity. 

4.4  Trickcr  Dat* 

Due  to  the  intrinsically  low  signal  sensitivity  of  the  on-board  receiver,  and 
to  tracker  hunting  errors,  large  fluctuations  occurred  in  the  instantaneous  position 
measurements  of  azimuth,  elevation,  and  slant  range.  These  data  were  refined 
after  each  flight  on  a  CDC  6600  computer.  The  data,  which  were  recorded  every 
0.  1  sec,  was  edited  by  passing  through  a  first  difference  software  routine.  Dif¬ 
ferences  between  adjacent  points  greater  than  two  standard  deviations  from  the 
average  deviation  were  flagged.  The  average  first  difference  is  defined  by: 


£  vy 

i=2 


n-1 


Where  n  is  the  number  of  points  and  X  ^  is  the  value  of  the  ith  point.  'The  standard 
deviation  is: 


1  2 


To  smooth  through  the  flagged  data,  a  search  was  made  to  find  four  statistically 
acceptable  points  on  either  side  of  the  flagged  values.  Once  found,  least  squares 
polynomials  were  generated  through  this  eight  point  array.  Root-mean-square  and 
Durbin-Watson  criteria  were  used  to  automatically  select  the  optimal  fit,  up  to 
fifth  degree.  Polynomial  interpolation  then  replaced  the  flagged  data. 

Smoothing  of  data  within  ±  2  o  was  performed  by  replacing  centered  data  with 
a  3  X  3  moving  average  of  neighboring  points.  The  average  is  defined  by: 

X  =  l/9X._2  +  2/9X._1  +  1/3X.  +  2/9X.+  1  +  l/9Xl+2  . 

Simultaneously,  a  second  arithmetic  running  average  was  performed  on  the 
resulting  data  to  create  10-sec  averages  spaced  at  2 -sec  intervals.  Furthermore, 
the  data  went  through  a  low -pass  digital  filter  having  a  -3  dB  band  edge  at  0.  01  Hz 
that  drops  to  -60  dB  at  0.  02  Hz. 
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The  smoothed  azimuth,  elevation,  and  slant  range  were  then  transposed  to  X, 
Y,  and  Z  Cartesian  and  geodetic  coordinates.  Using  data  time  spacings  of  30  sec, 

north,  east,  vertical,  horizontal  balloon  speeds,  and  horizontal  wind  shears  were 
calculated. 

Sample  results  of  the  balloon  position,  velocity,  and  horizontal  wind  shear 

calculations  are  presented  in  Figure  7.  The  data  are  interpolated  for  even  spac- 

£ 

ings  of  time  every  2.  0  second. 

5.  PROGRAM  DESCRIPTION 


The  general  resear cr.  objective  of  this  program  is  to  compare  the  measure¬ 
ments  of  atmospheric  turbulence  of  several  different  experimental  designs.  The 
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primary  instruments  consisted  of  the  MIT  Lincoln  Laboratory  440  MHz  turbu- 
lence  scatter  radar  and  an  Air  Force  Weapons  Laboratory  Stellar  Scintillometer.  ^ 
In  addition  to  the  primary  instrumentation,  radiosonde,  thermosonde,  and  air¬ 
craft  instrumentation  were  used  to  provide  a  broader  comparison  base,  AFGL 

operated  the  sondes  and  scintillometer,  while  personnel  from  the  Northeast  Radio 
Observatory  Corporation  (NEROC)  operated  the  radar.  Aircraft  m  easurements 
were  provided  by  Airborne  Research  Associates,  Bedford,  Mass.  ,  under  contract 
to  AFGL. 

From  10  Feb  1081  until  8  Aug  1081,  over  30  radiosondes  were  flown.  Twenty- 
nine  of  these  are  listed  in  Table  1.  The  preliminary  qualification  flights,  up  until 
21  May  1981,  were  Launched  from  AFGL.  The  remainder  of  the  experiments  (22) 
were  launched  near  the  NEROC  Haystack  Observatory  at  Millstone  Hill,  Westford, 
Mass.  ,  during  the  month  of  July. 

Scheduling  conflicts,  time  of  day,  and  equipment  availability  dictated,  to  some 
extent,  the  particular  experiment  on  a  given  day.  The  scintillometer  could  be 
operated  only  on  cloudless  nights.  The  aircraft  required  co-ordination  with  the 
FAA,  and  it  was  further  constrained  by  bad  weather  and  other  scheduling  difficulties. 
The  radar  was  shared  by  other  groups  at  Millstone,  and  it  was  available  for  these 
measurements  only  on  those  dates  pre-arranged  weeks  in  advance.  The  radiosonde 

■'As  of  this  printing,  a  detailed  description  of  the  positional  data  was  being 
prepared  for  publication  as  an  AFGL  technical  memorandum.  The  title  is: 

Analysis  of  Balloon  Trajectory  for  the  Calculation  of  Wind  Velocities  and  Shears, 
by  Sullivan,  B. ,  et.  al. 

12.  Watkins,  B.J.  (1081)  Millstone  Hill  Doppler  Radar  Observations  of  Winds  and 

Turbulence,  Middle  Atmospheric  Program, "Handbook  for  NR]?,  Extended 
abstracts  from  Int.  Sym.  on  Middle  Atm.  Dyn.  2:30-36. 

13.  Ochs,  G.R.,  Lawrence,  R.S.,  and  Wang,  T.  (197  6)  Stellar  scintillation 

measurement  of  the  vertical  profile  of  refractive— index  turbulence  in  the 
atmosphere,  Proc.  Spc.  of  Phwto-optical  Inst.  Eng.  75:48-54. 
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Fig-are  7.  Sample  Output  of  Reduced  Tracker  Data.  Edited  and  filtered 
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Table  1.  Optical  C^,  (Turbulence)  Measurement  Program  (Continued) 


and  thermosonde  had  no  constraints  other  than  availability  of  tracking  personnel. 

As  it  turned  out,  the  scintillometer  and  radar  operated  together  only  on  one  even¬ 
ing,  31  July  1981.  Good  quality  experiments  using  three  techniques,  radar, 
thermosonde,  and  aircraft  were  accomplished  twice.  Good  quality  radar  and  radio¬ 
sonde  measurements  were  completed  on  fifteen  occasions.  Good  quality  radiosonde 
and  thermosonde  experiments  were  operated  on  eighteen  occasions.  Radar,  radio¬ 
sonde,  and  thermosonde  operated  successfully  on  eleven  occasions. 


6.  RESULTS 

Results  of  typical  experimental  runs  are  presented  in  the  following  plots. 
Figures  8  and  9  show  temperature  and  relative  humidity  profiles  obtained  31  July. 
A  small  temperature  inversion  is  found  at  3  km  that  roughly  corresponds  with  a 
high  negative  gradient  in  relative  humidity.  A  more  pronounced  temperature 
inversion  is  found  at  12  km,  defining  the  tropopause. 


Q/20  MET  OPT ft 


Figure  8.  Typical  Temperature  Profile  Plot,  Processed  result 
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Figure  9.  Typical  Relative  Humidity  Profile  Plot.  Processed 
result 


Figures  10  and  11  show  the  east  and  north  balloon  velocity  component  profiles. 
The  maximum  horizontal  wind  occurs  at  10.7  km  at  35  m/sec.  The  detailed  wind 
structure  is  clearly  evident. 

Figure  12  shows  the  horizontal  vector  wind  shear,  S^,  profile  to  30  km. 

is  defined  as  the  resultant  of  the  east  and  north  horizontal  component  shears,  or 
2  2  2 

SH  =  SN  +  S£.  Figure  13  is  an  expanded  view  of  Figure  12  in  the  altitude  range 
from  20  km  to  25  km  and  shows  the  mean  shear  at  about  0.02  sec  .  This  plot 
shows  structuring  of  the  wind  shear  down  to  the  spatial  filter  wavelength  of  about 

120  m.  The  power  spectral  density  of  the  winds  are  plotted  in  Figure  14.  The 

14 

slope  is  about  -2.  5,  which  compares  well  with  other  measurements. 


14.  Endllch,  R.M.,  and  Singleton,  R.C.  (1969)  Spectral  analysis  of  detailed 
vertical  wind  speeds  profiles,  J.  Atmos.  Scl.  20:1030-1041. 
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KH! 


Figure  ll.  Typical  Balloon  Velocity  Analysis 
North  Profile  Plot.  Edited  and  filtered 


OXOQFO  FINAL  TWCAT  ANALYSIS 


$H, VERTICAL  SHEAR  Or  HORIZONTAL  WIND 

Figure  12.  Typical  Horizontal  Vector  Shear  Analysis. 
Profile  plot.  Shear  is  calculated  over  a  30-sec  data 
spacing 

QXOOK3  FINAL  TRAOAT  ANALYSIS 


SH, VERTICAL  SHEAR  OF  HORIZONTAL  WINO 

Figure  13.  Expanded  View  of  Horizontal  Vector  Shear. 
The  shear  altitude  resolution  is  about  12  0  m 
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Figure  14.  Power  Spectral  Density  Analysis  of  Horizontal 
Vector  Shear.  The  slope  is  about  -2.  5 


Figure  15  shows  "RAW"  thermosonde  C^,  values  as  a  function  of  altitude.  The 

mean  tends  to  decrease  logarithmically  with  altitude.  Figure  16  is  an  expanded 

view  of  Figure  15  in  the  altitude  range  from  10  km  to  15  km.  A  large-scale  max- 

mum  shows  up  at  11.  5  km  while  the  fluctuations  in  C.'  -  appear  to  be  due  to  the  small 
2 

sampling  interval.  C.,  smoothed  to  120  m  is  plotted  in  Figure  17  for  comparison. 

2 

Figure  18  plots  the  cumulative  frequency  distribution  of  C^,.  The  data  follows  a 
straight  line,  except  at  the  endpoints,  which  confirms  the  log  normal  distribution 
of  C2n. 

Tlie  purpose  of  the  optical  turbulence  program  was  to  compare  various  tech- 
2 

niques  of  observing  C^j.  Sonde  measurements  were  Important  in  the  low  altitude 

comparison  between  radar,  scintillometer,  and  thermosonde.  The  radar-received 

2 

power,  scattered  by  refractive  index  fluctuations,  can  be  related  to  a  Cjj.  This 
is  caused  by  turbulent  fluctuations  in  both  temperature  and  water  vapor  concen¬ 
trations.  As  the  comparisons  are  based  on  optical  wavelength  C^,  the  contribution 
of  water  vapor  must  be  removed  from  the  radar  C^. 
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1  15 

The  fluctuations  of  the  refractive  index  across  a  layer  is  given  '  as: 


M  «  -77.6  X  10 


“  L  WB 


where 


P  *  pressure  (mb) 

T  *  temperature  (Kelvin) 

_2 

g  ■=  gravity  (m  *  s  ) 
q  =  specific  humidity  (gm/gm) 
u-2  ■  squared  Brunt-Vaisala  frequency  (s  2)  '  (g/6)®1 
q1  =  vertical  gradient  of  specific  humidity  (m  ) 

6'  =  vertical  gradient  of  potential  temperature  (Kelvin/m) 

2  2 
The  radar  C^,  can  be  corrected  to  dry  conditions,  (dry)  by  evaluating  the 

bracketed  terms  in  Eq.  (8)  on  the  basis  of  the  sonde  measurements  of  P,  T.  and  q. 


C^j(dry)  *  C^j(radar) 


1  +  15(0oo^  -H^OOJL^I5 


The  comparison  between  the  CN  (dry)  derived  from  radar  measurements  and  the 

optical  scintillometer  is  shown  in  Figure  19.  The  optical  scintillometer  obtains 

-lata  over  very  large  weighting  functions.  The  scintillometer  i/e  naif-width  varies 

from  1.4  km  at  an  altitude  of  2.2  to  4.  9  km  at  an  altitude  of  14  km.  The  weighting 

functions  are  shown  in  Figure  20.  The  comparison  between  radar  and  scintillom- 

2 

eter  shown  in  Figure  19  is  obtained  by  passing  the  (dry)  through  the  six  weight¬ 
ing  functions . 

A  comparison  of  thermosonde  data  and  the  scintillometer  is  shown  in  Fig¬ 
ure  2  1.  Here  the  thermosonde  data  shown  in  Figure  15  is  passed  through  the 
scintillometer  weighting  functions.  This  agreement  is  good.  The  large  smoothing 
of  the  weighting  functions  hides  whatever  variations  may  exist  due  to  the  space  and 
time  differences  in  the  respective  measurements.  The  scintillometer  viewed  the 
sky  vertically  while  the  thermosonde  followed  a  slant  path. 

Thermosonde  and  radar  comparison  can  be  made  over  very  nearly  the  same 
volume  and  space  by  the  fact  that  the  radar  can  be  pointed  in  the  direction  of  the 
sonde.  The  radar  (dry)  is  obtained  over  one -minute  periods  and  averaged  over 


15.  VanZandt.  T.  E.  ,  Gage,  K.S.,  and  Warnock,  J.M.  (1981)  An  Improved  Model 
for  the  Calculation  of  Profiles  of  and  e.  in  the  Free  Atmosphere  from 
Background  Profiles  of  Wind,  Temperature,  and  Humidity,  2  0th  conf.  on 
Radar  Meteorology,  (3!)  Nov  — 3~Bec  198 1). 
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five  minutes.  A  composite  of  the  radar  results  is  made  to  match  the  sonde  altitude 

2 

and  times.  The  comparison  of  radar  and  thermosonde  values  is  presented  in 
Figure  22.  The  thermosonde  data  has  been  averaged  over  0.  50-km  Intervals. 

The  comparisons  shown  Indicate  the  thermosonde  radar  and  scintillometer 
are  in  quite  good  agreement.  The  data  shown  represents  the  results  from  just 
one  day's  measurements.  Detailed  comparisons  involving  all  the  data  sets  will 
be  analyzed  and  presented  in  separate  acientific  reports. 
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Figure  15.  Raw  Thermosonde  Analysis  -  Profile  Plot.  Low  and 
high  gain  measurements  have  been  combined  into  a  single  best  meas¬ 
urement.  The  thermosqnde  measures  Cyj  from  ground  level  to  30  km 
CL  is  calculated  from  and  the  radiosonde  measurements  of  pres¬ 
sure  and  temperature 
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ALT  IT 


ALTITUDE  (km) 


Figure  10.  Comparison  of  Radar  and  Scintillometer  C 
Profiles.  The  path  averaged  radar  and  scintillometer  pro 
files  show  good  agreement.  The  points  depict  radar 
Cjj  measurements.  The  solid  line  through  the  radar  points 
represent  5 -min  radar  averages.  The  seven  scintillometer 
measurements  are  shown  as  *  with  straight  lines  drawn 
between  them.  The  dashed  curve  represents  the  radar 
averaged  C?j  corrected  for  moisture  and  passed  through  the 
scintillometer  weighting  functions  of  Figure  19 
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Figure  21.  Comoarlson  of  Thermosonde  Cjg  and 
Scintillometer  Cfj  and  Scintillometer  Cjj  Profiles. 
The  Thermosonde  "raw"  Cfj  data  was  passed 
through  the  scintillometer  weighting  functions  of 
Figure  18.  Although  the  scintillometer  viewed  the 
sky  vertically  while  the  thermosonde  followed  a 
slant  path,  the  agreement  is  good 
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ALTITUDE  (km) 


igure  22.  Comparison  of  Ther™ *  The  R&S  dite^s 
pace  coordinates  of  the  Radar  and  Thermosonde, 


TIME  (balloon) 
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